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ABSTRACT: The symmetry constraints imposing on the quantum states 
of a dot with 13 electrons has been investigated. Based on this study, the 
favorable structures (FSs) of each state has been identified. Numerical cal- 
culations have been performed to inspect the role played by the FSs. It was 
found that, if a first-state has a remarkably competitive FS, this FS would 
be pursued and the state would be crystal-like and have a specific core-ring 
structure associated with the FS. The magic numbers are found to be closely 
related to the FSs. 

PACS (numbers): 73.61.-r 

1, INTRODUCTION 

Modern experimental techniques, e.g., by using electrostatic gates 
and by etching, allow a certain number of electrons to be confined in semi- 
conductor heterostructures. 1-6 Such many-electron systems have much in 
common with atoms, yet they are man-made structures and are usually called 
" quantum dot " . The number of electrons contained in a dot ranges from a 
few to a few thousands, they are confined in a domain one hundred or more 
times larger than the atoms. Thus, in addition to atoms, nuclei, - • • that exist 
in nature, quantum dots as a new kind of system will definitely contain new 
and rich physics, and therefore they attract certainly the interest of academic 
research. 

On the other hand, the properties of the dots can be changed in 
a controlled way, e.g., by changing the gate voltage or by applying an ad- 
justable magnetic field, etc. Therefore, these systems have a great potential 
in application. Due to this fact, the investigation of quantum dots is a hot 
topic in recent years 1-6 . In the experimental aspect, progress has been 
made to reveal different kinds of physical property. A crucial point is to 
clarify the electronic structures. An important step along this line is the 
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first observation of the Coulomb blockade spectra via the measurement of 
conductance as a function of gate voltage 7 , where very clear level structure 
has been demonstrated. Afterwards, a substantial amount of information is 
drawn from conductance measurement. The measurement of the difference 
in chemical potential exhibits also clear shell structures 8 . The excitation 
of electron can be probed by far-infrared and capacitance spectroscopy 9,10 
With further progress in experimental techniques, the dots will definitely be 
understood better and better, and they will serve as a rich source of infor- 
mation on many-body physics in the coming years. 

In the theoretical aspects, detailed information on electronic struc- 
tures has been obtained for the systems with a smaller N (say, N<10) 2,4 ' 6 . 
When N is small, the effect of symmetry was found to be very important, 
e.g., the magic angular momenta of few-electron dots originate from the 
constraint of symmetry 11-13 . When N is larger (say, N >10), the effect of 
symmetry is scarcely studied. The systems with a larger N are themselves 
very attractive, because they might possess both the features of few-body 
and many-body systems. Thus the understanding of these systems might 
serve as a bridge to connect few-body and many-body physics. In a previ- 
ous paper, the electronic structures of a dot with nine electrons have been 
studied 13 . The present paper is a continuation of the previous one, and is 
dedicated to the study of the dot with N=13 and with the spins polarized. 
The choice of thirteen is rather arbitrary, just because it is explicitly larger 
than the systems with N<10 which have already been extensively studied, 
and because it is not very large so that accurate numerical calculations (in 
the qualitative sense) and detailed analysis can still be performed. From a 
previous study by a number of authors 13-18 , it is believed that a general 
picture of dots would consist of a core surrounding by a ring. It would be 
interesting to see, when N is larger, how the details of the core-ring structure 
would be and how these structures are affected by symmetry . Such a study 
would exhibit further insight of many-body physics. 

In what follows, the 13-body Schrodinger equation is solved via an 
exact diagonalization of the Hamiltonian, the accuracy has been evaluated. 
The underlying dynamical and symmetry background has been studied. Fa- 
vorable structures for each state have been suggested based on symmetry 
consideration. The eigenwavefunctions have been analyzed in detail to ex- 
hibit how the electronic structures are affected by symmetry. The appearance 
of magic numbers is discussed. 
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2, HAMILTONIAN AND THE APPROACH 

Let the electrons be fully polarized (therefore the spin-part can be 
neglected and the spatial wave functions are totally antisymmetric) , and con- 
fined in a 2-dimensional plane by a parabolic confinement. The Hamiltonian 

reads 

H = T + U (1.1) 
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where m* is the effective electron mass, e r is the dielectric constant, and 
fta> measures the strength of the parabolic confinement ( uj arises mainly 
from a magnetic field B. This field leads also to a term linearly proportional 
to B . This term has been neglected because it does not at all affect the 
eigenwavefunctions, and therefore not affect the electronic structures). 

In order to diagonalize the Hamiltonian, a set of orthonormalized 
single-particle harmonic oscillation (h.o.)states <p m k are introduced. Here, 
<p m k is an eigenstate of a pure h.o. Hamiltonian 

h = -2&V 2 +|m*fi2 r 2 (2) 

where Qq is an adjustable parameter in general not equal to u , This 
eigenstate has eigenenergy (m + k + l)hQo and angular momentum (m — k)h 
. From them the many-body basis functions (BFs) 
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with a given total orbital angular momentum L = Ei( m i — h) are com- 
posed. 

From the BFs, the eigenstates of the dot are expanded as 

* = E^ Q (4) 
where the coefficients C a can be obtained via a procedure of diago- 
nalization. The ip a are arranged in such a sequence that < ip a \H\ifj a > < 
< ip a+ i\H\ ip a+ i > . Evidently, in such a sequence, the ifj a with a smaller 
a is more important to the low-lying states, while those with a very large a 
can be neglected. The H will be diagonalized step by step. In the first step, 
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H is diagonalized in a smaller space with N a BFs ( ipi to ipN a ) ■ Then, H 
is diagonalized again in a larger space with N& BFs ( from tpi to ipN b , and 
N b is considerably larger than N a ) . This process repeats again and again 
until a satisfactory convergency of the lower eigenenergies is achieved. In the 
first step, all the ip a for the diagonalization is limited to the lowest Landau 
levels (LLL), i.e., all the <fi mi ki contained in ip a have hi = 0. However, step 
by step, BFs of higher Landau levels will mixed in. In order to speed up 
the convergency, the flo in eq. (2) is considered as a variational parameter to 
optimize the lower eigenenergies emerged from the diagonalization. 

In the following calculation, we have m*=0.067m e , hu)o=3meV, 
e r =12A (for a GaAs dot). To show the convergency, as an example, the low- 
est eigenenergies with L=82 are obtained as 436.895, 436.806 and 436.760meV 
when the number of BFs are 6000, 9000, and 12000, respectively. One can see 
that the convergency is not very good. However, the densities calculated be- 
low by using 6000, 9000 and 12000 BFs are indistinguishable (e.g. in Fig.l). 
Since we are mainly interested in the qualitative aspect, the accuracy that 
we have achieved is sufficient. 

After the diagonalization the eigenstates are obtained. The series 
of states having the same L is labeled as (L)j. The i = 1 state (the lowest 
of the L-series) is called a first-state. 

The eigenwavefunctions of a 13-electron system are complicated. In 
order to extract informations from them the following physical quantities are 
defined and calculated. They are the one-body density 

pi(ri) = / \^ L \ 2 dr 2 dr 3 ■ ■ ■ dr l3 , (5a) 

the two-body density 

p 2 (ri,r 2 ) = / l^il 2 dr 3 (ir4 • • • dr 13 , (5b) 
and the three-body density 

P3(ri,r 2 ,r 3 ) = / \^ L \ 2 dr 4 dr 5 ■ ■ ■ dr 13 , (5c) 

It was found that in many cases the pi(r) has an outer peak and an 
inner peak with a minimum lying in between (at r = a). In this case we can 
define an outer region (r > a) and an inner region (r < a) . Accordingly, 
we can define the average number of particles N out and N in contained in the 
outer and inner regions, respectively, as 

N out = Nf~p 1 (r 1 )dr 1 (6a) 
N in = NfZp 1 (r 1 )dr 1 (6b) 

For example, the (88)i state has a =367.81, N out = 9.97 and N in =3.03. 
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Once the border a is defined, we can define the angular momenta l ou t 
and the moments of inertia I out contributed by the outer region, respectively 

, as 

lout = N f™ dn / ** L h*Ldr 2 ■ ■ ■ dr 13 (7) 

and 

I out = M/ a °°Pi(ri)r^r 1 (8) 
where M = Nm* is the total mass. Similarly, the l in and I in contributed 
from the inner region can also be defined. Although these quantities are not 
good quantum numbers, they can help us to understand better the physics 
as shown later. 

3, DYNAMICAL AND SYMMETRY BACKGROUND 

Quantum mechanic systems are affected by both dynamical reasons 
and symmetry consideration . The following points are noticeable: 
(i) Core-ring structures. 
The spatial wave functions of low-lying states are mainly distributed in an 
area where the total potential energy U (eq.(1.3)) is lower. In particular, they 
would like to be distributed surrounding the (local) minima of U. In order to 
find out the (local) minima, let N in electrons be contained inside to form a 
core, and N out electrons be contained outside to form a ring, N in +N ou4 =N. 
When the relative locations of the electrons are appropriately adjusted (e.g., 
they form two homocentric regular polygons with or without an electron at 
the center) U will be optimized and arrives at its (local) minimum U op t, the 
associated configuration is called an N in — ~N out core-ring configuration. In 
this configuration, let the ratio of the radii of the outer polygon and the 
inner polygon be denoted as G opt . U opt and G opt are given in Table 1. 

Table 1, The optimal values U op t and the associated G op t of the (lo- 
cal) minima of U , each is associated with a N in — N out core-ring configuration. 



Nj re N otl £ 


1-12 


2-11 


3-10 


4-9 


5-8 


6-7 


7-6 


8-5 


U opt (meV) 


281.41 


278.35 


274.83 


274.22 


274.93 


276.31 


275.73 


278.76 


G op t 




3.73 


2.83 


2.41 


2.16 


2.02 


1.77 


1.70 



Evidently, a too small or too large N out (say, N out < 7 or N out > 11) is not 
advantageous to binding. Furthermore, the outer polygon should be neither 
too close to nor too far away from the core. 

In what follows, when the wave function of a state is distributed 
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surrounding a Nj n — N OU ( core-ring configuration, then the state is said to have 
a Nj„— N OU i structure. If the configuration has an electron at the center, then 
the structure is further denoted as (N in ) c — N out , otherwise as (Nj„) fe — N out 
. The subscript h implies a hollow structure. 

It is shown in the table that the U opt of a number of configurations 
are quite close to each other. At a first glance, one might expect that a 
strong mixing of geometric configurations would occur and would spoil the 
crystal-like picture. However, this is not true mainly due to the quantum 
constraints as we shall see later. 

(ii) Uniform rotation . 

Let us consider first a classical model system of two rotating ho- 
mocentric rings, the outer ring has (b < r < a), while the inner ring 
has (d < r < c , and c< b). Let the angular momentum, the mo- 
ment of inertia and the angular velocity of the outer (inner) ring be l out , 
I out and ijj out ( l in , I in and u>i n ), respectively. The total angular momen- 
tum L = 1^ + lin = lout^out + hnUin , and the total rotation energy 
T = \(I ut^lut + hn^fn) ■ Now, let us ask how the u out and u in would 
be chosen so that T is minimized under the condition that L is conserved? 
The answer is simply u ou t = Win — L/(I out + I in ) = L/I. This fact implies 
that if the two rings are rotating with the same angular velocity, the rota- 
tion energy can be reduced. Although this point is viewed from classical 
mechanics, however the first-states of a quantum mechanic system would do 
its best to lower the energy, thus they would pursue a uniform rotation, i.e., 

Wout ~ Wi n ■ 

From the point of view of quantum mechanics, the low-lying states 
are mainly dominated by the BFs belonging to the LLL. In these BFs, all the 
single-particle state (f> m k have k = and angular momentum I = m — k = m. 
For each 4> mk , the angular velocity can be defined as uo =< I > / (m* < r 2 >), 
which is proportional to ^ if k — 0. Evidently, uj is close to a constant 
unless I is very small. Thus , for the BFs of the LLL, all the electrons rotate 
with similar angular velocities, and we have the uniform rotation u ou t ~ Wi n 

(iii) Symmetry constraints and the favorable structures. 

It has been found that inherent nodal surfaces are imposed in wave func- 
tions by symmetry, thereby the structures of quantum states are seriously 
affected. 12,19 ^ 21 In the case of 2-dimensional polarized quantum dots , it was 
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found that a wave function would be zero when the electrons locate at the 
vertexes of a regular N-side polygon if L 7^N(j + ) , where j is an 

integer 6,11,22 . This constraint can be generalized to the core-ring structures. 
Let the ring has an angular momentum l out , while the core has li n . When 
the outer particles locate at the vertexes of a N out -side polygon, and the in- 
ner particles locate at the vertexes of a N D -side polygon (N G =Nj n or N in — 1, 
in the latter case an electron would stay at the center) , then it is straight 
forward to prove that the wave function would be zero if 

lout ^out(j2 + (1 + {-l) N ^)/A) (9a) 

or 

kn ^KUi + (! + (-1)^/4) (96) 

where ji and j'2 are integers. In other words, the above configuration 
would be prohibited if l out ( l in ) does not relate to ~N out (N Q ) in the above 
way. Thus, a (N in )—N out structure would be pursued by a first-state only if 
the L can be divided as a sum of li n and l out so that the requirements (9a) 
and (9b) are fulfilled. If this happens, the (Nj n )— N out structure is called 
a candidate of favorable structure ( CFS) of the state. Incidentally, for an 
eigenstate, both the l out and U n are not good quantum numbers, they appear 
as the angular momenta of the main component of eigenwavefunctions. 

Usually each state may have a number of CFS, some of them 
are not competitive due to having a too small or too large N out , they can be 
neglected. In what follows, among the CFS of a state, if some of them have 
N ou < > 8, then those with N out < 7 are neglected; if all the CFS has N out < 7, 
then all of them would be neglected except the one with the largest N out ; all 
the CFS with N out = 12 are neglected without exception. After the neglect, 
the remaining CFS are call the favorable structures (FSs), they are listed 
in Table 2. E.g., the L = 86 state has four CFS , namely the (5) c -8, (9) c -4, 
(ll) c -2, and (12) c -l . Among them only the (5) c -8 as a FS is listed in Table 
2. When a state has more than three FS, only the most competitive three 
are listed. 

(iv) Excitation of the core 

This paper concerns only the low-lying states with L > N(N — 1) /2 ( 
or the filling factor v < 1 ), they contain mainly the BFs belonging to the 
LLL. In these BFs, the angular momenta of any pair of electrons can not 
be the same due to the Pauli Principle. Therefore, they can be denoted as 
ipa = {lik ■ ■ - In} with U < 
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Since we have 

< Mr 2 \(t>io >= (I + (10) 
the spatial distribution of the wave function ro depends on I. The 
smaller the /, the smaller the size. Thus, for the ip a belonging to the LLL, 
U n is just equal to l\ + 12 + • • • + Zjv in ■ If h — 0, there must be an electron 
staying at the center, because the 0oo wave function is distributed closely 
surrounding the center. Thus, a (Nj n ) c -N out structure must be contributed 
by the ip a with li — , while a (N in ) h -N out structure is contributed by those 
with li > . When all the k of the inner electrons satisfies k + 1 = , 
the inner electrons are said to be compactly aligned. Meanwhile, Li n would 
arrive at its lower bound (h n )b =N in (Nj„ — 1) /2 , if l\ = 0. In this case, we 
say that the core is inert (not excited). Otherwise, we have l in > (h n )b, and 
we say that the core is excited. Evidently, all the hollow states must have 
li > 0, thus they have an excited core. 

When L < 90, core excitation is not possible (unless the electrons 
jump to higher LLL), therefore the first-states would have a (Nj„) c — N out 
structure with the core inert. However, when L > 90, core excitation might 
occur. It implies two cases: (a) The inner electrons have their l { remaining 
to be compactly aligned but with l ± — k , and therefore have a (N in )/j-N out 
hollow structure . (b) The U of the core are no longer aligned compactly, 
e.g., li — while l 2 = 2, etc. . 

It was found that, when L is not large (say, L < 101 ), the first-states 
have either an inert core or an excited compact core with l ± — 1 , as shown 
in Table 2 . However, when L is large, higher core excitation with li > 1 will 
emerge as shown later. 

Incidentally, due to eq.(10), the compact alignment of the angular 
momenta also implies a compact alignment of radial positions. Thus, in 
the core-ring structures, the groups of inner and outer electrons may each 
compactly aligned. The associated BF can be simply denoted as {li — 
lN in ,lN in +i ~ In} (e.g., {1,2,3, 6,7,- • -,15} ={1-3,6-15}.). This is called a 
two-bunched BF by Ruan 23 (The one-bunched BF {l± — is a special case 
of two-bunched BF with l^ in = £jv in +i — 1 )• It is straight forward to prove 
that, for a CFS of a v < 1 state, among the BFs of the CFS, one and only 
one of them is a two-bunched BF belonging to the LLL. Thus, a simple way 
to find out the CFS is to look for the two-bunched BFs of a state. 
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(v) Particle separation 

It is noted that the U in eq.(1.3) can be exactly rewritten as 

U = \Mu 2 R 2 c + E Uefffak) (11.1) 

j<k 

Where R c is the radial distance of the cm., and M e //( r jfc) is the effective 
pairwise interaction 

There is a minimum in u e ff located at rj k = r u = ( 47r£ e (m*u ) 2 )^ 3 • 
Evidently, if each electron separates from all its adjacent electrons by this 
distance, the potential energy can be minimized. Therefore, in low-lying 
states, adjacent electrons would roughly keep the separation r u . With the 

o 

above parameters, r u = 576.3 A .The r u is a basic measure and is useful 
for the understanding of electronic correlation and the size of the system. 
Obviously, for a (N in )— ~N out structure, the distance between the ring and the 
core depends closely on r u . 

(vi) Core-ring separation 

It is recalled that, in order to minimize the potential energy, the ring 
should separate from the core by an appropriate distance. In this subsection, 
we shall evaluate the core-ring separation by using the approximation of 
uniform rotation. 

For a given CFS with the given N out and l out , let us define a quantity 

n / lout/tin / 1 H\ 

y - V N out /N in y iZ ) 

On the other hand, we have 

lout — m *N out < T 2 >ri ng (13) 

(this equation is the same as eq.(8)), and a similar equation for I in . Thus 
we have 

9— \j l ou //l" n ' (< r<1 >ring / < T 2 > core ) = '• (< T 2 > r ing / < T 2 > C ore) 

It is believed that the uniform rotation is a good approximation for the 
first-states, because they should do their best to lower the energy (this is a 
point remain to be checked). Under this approximation 
9k (< r 2 > ring I < r 2 > core fl 2 (15) 

The optimal value of the right hand side of eq.(15) has been denoted as 
G op t given in Table 1. Thus, if a FS has its 9 (evaluated from the definition 
eq.(12)) close to G opt , then the core-ring separation is appropriate and the 
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FS is advantageous to binding and therefore competitive. Otherwise, it is 
not. 

The g /Gop t of the FS are also listed in Table 2, many of them are found 
to be very close to one. E.g., the FS of the (86)1 is a (5) c — 8 structure with g 
=2.18, the associated G opt is 2.16 (cf. Table 1) , thus they are close to each 
other . 

The above points are important to the following discussion. When a 
first-state has a FS which is superior than the other FSs ( or the state has only 
one FS), the FS is expected to be dominant. In this case the state would have 
a clear geometric feature arising from the N in —N out structure of the FS, and 
appear to be crystal-like. However, when a first-state has a few nearly equally 
competitive FSs, its structure can not be uniquely predicted. Nevertheless, 
the Table 2 is a key to understand the electronic structures. 

Table 2, Characters of the first-state from symmetry consideration 
and from our calculations. 
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Favorable structures and 




Quantities evaluated from 




their related features 




the pi 


of the first-states 


L 


FS 


h 


kn 


9 


9/ G opt 


a 


N m 


kn 


7 


81 


(10) c -3 





45 


1.6330 












82 


(9) c -4 





36 


1.6956 












83 


(8) c -5 





28 


1.7728 


1.04 










84 


(7)c-6 





21 


1.8708 


1.05 


2.700 


6.55 


21.30 


0.89 


85 


(6) c -7 





15 


2.0000 


0.99 


2.550 


5.76 


16.45 


0.91 


86 


(5) c -8 





10 


2.1794 


1.01 


2.288 


4.65 


10.52 


0.93 


87 


(4) c -9 





6 


2.4495 


1.02 


2.100 


3.85 


7.21 


0.95 


88 


(3) c -10 





3 


2.9155 


1.03 


1.889 


3.03 


4.44 


0.96 


89 


(2) c -ll 





1 


4.0000 


1.07 


1.555 


2.00 


1.90 


0.98 


90 


(l) c -12 












1.120 


0.97 


0.42 


1.094 


92 


(6)c-7 





15 


2.0976 


1.04 


2.700 


6.26 


19.70 


0.94 


93 


(8) c -5 





28 


1.9272 


1.13 


2.414 


4.99 


12.43 


0.964 


94 


(5) c -8 





10 


2.2913 


1.06 


2.377 


5.01 


12.51 


0.967 


95 


(9) h -4 





45 


1.5811 




2.205 


4.06 


8.30 


0.967 


96 


(4) c -9 





6 


2.5820 


1.07 


2.181 


3.89 


7.38 


0.998 


97 


(7)/»-6 


1 


28 


1.6956 


0.96 


1.942 


3.04 


4.80 


0.967 


98 


(3) c -10 





3 


3.0822 


1.08 


1.926 


2.91 


4.19 


1.012 


99 


(5) fc -8 


1 


15 


1.8708 


0.86 


1.519 


1.74 


1.86 


0.904 


100 


(2) c -ll 





1 


4.2426 


1.14 


1.611 


1.90 


1.98 


0.981 


100 


(4)^-9 


1 


10 


2.0000 


0.83 










101 


(3) h -10 


1 


6 


2.1794 


0.77 


1.936 


2.82 


6.54 


0.813 



4, EIGENENERGIES 

After performing the diagonalization, eigenenergies and eigenstates 
are obtained. Let E((L)i) be the energy of the (L), state. It is noted 
that , for a first-state, if the Coulomb repulsion among the electrons are 
removed, all the electrons would fall in the LLL with the energy (L + N)hu 
. For this reason, let us define e(L) = E((L)i) — (L + N)hu .This quantity 
is a measure of the Coulomb repulsion in the first-states, which is plotted 
in Fig. 2 in accord with L. When L increases, the size of the system will 
increase a little , the Coulomb repulsion will thereby decrease. Thus, s(L) 
decreases monotonously with L as shown in the figure. However, there are 
four platforms. We shall return to this point later. 

5, ELECTRONIC STRUCTURES (78 < L < 90) 
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In what follows mainly the results of the first-states are given. We 

use dM = \J — ir^- = 194. 71 A as the unit of length. The optimal separation 
r u = 2.96a M . 

Let us begin from the state with the filling factor v — 1 , namely 
the (78)i state. This state has only one BF {0,1,2,- • -12} (for short, {0-12}) 
belonging to the LLL, this BF has a weight 85.5%. In this BF, the electrons 
are roughly uniformly distributed inside a circle as shown in Fig. 3a. It is 
noted that a clear geometric structure arises from the coherent mixing of 
BFs. Due to the lack of mixing, the (78) i can not have a clear geometric 
structure, therefore it is liquid-like as shown in Fig. 4a. 

On the other hand, for the number N together with two arbitrary 
integers n (<N), and f, there is an identity 

N(N-l) + -, N = n(n+2j'-l) + (N-n)(N+n+2f -1) ^ 

Let the left hand side be equal to L, and the two terms at the right be 
equal to k n and l out . Then this identity is associated with a division of L. 
When f — , the left hand side of (16) is equal to 78 . It is easy to see 
that the pair N G = n — 1 and U n meet the requirement of eq.(9b), while the 
pair N out =N — n and l out meet the requirement of eq.(9a). Thus, eq.(16) 
implies that all the (N G + l) c -N 0Ui structures with N G = to 12 are the CFS 
of the L = 78 states. Therefore the wave function of the (78)ican get access 
to all the symmetric configurations 12,13 , and thus is nodeless (except a pair 
of electrons overlap with each other). Accordingly, the energy of this state 
is lower. 

For the (79)i state, there is also only one BF {0-11,13} belonging to 
the LLL. Thus this state is also liquid-like as shown in Fig. 4b. However, on 
the contrary with the (78) i , all the (Nj n )-N OU ( structures are not the CFS of 
the (79)i, except the (12) c -l which is very poor in binding. Thus the energy 
of this state is much higher. Owing to the (78)i is lower while the (79)i is 
higher, the difference leads to a platform appearing in Fig. 2 between L = 78 
and 79. 

Ranging from (79)i to (90)i, all these states have only one FS, thus 
their structures can be well predicted. The N out of their FS (cf. Table 2) 
increases from 1 to 12, this leads to a regular variation of their electronic 
structure. When N ou< is small (say, N out < 5 ), the outward electrons are 
found to be very close to the core. As a result, their ring-core-structures are 
ambiguous as shown in Fig. 3b and 4b, where the patterns are representative 
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for the (79)i to (83) 1 states. In these states the FS itself is not competitive. 
This fact would lead to a stronger mixing of structures, and therefore they 
are liquid-like. 

Even in the liquid-like states, electronic correlation can still be 
viewed via the three-body densities as shown in Fig. 5a and 5b, they are 
representative. Fig.5a for the (81)i exhibits that the three outward electrons 
(two are labelled by white spots and one by a double-peak, which implies an 
oscillation around an equilibrium position) are very close to the core. This 
fact supports the presumption that the FS, namely the (10) c -3 structure (cf. 
Table 2), is pursued by the state . Although the U of the (10) c -3 is higher, 
however no other better symmetric configurations are allowed by symme- 
try. Consequently, the component of the (10) c -3 is still relatively important 
Since the outward electrons are so close, the core is strongly deformed. 
There are three peaks at the outer ridge of the core, it implies that three 
inward electrons form a regular triangle close to the border. Fig. 5b for the 
(82) 1 exhibits that the four outward electrons are also very close to the core. 
This fact supports again that the FS is pursued. The core is also strongly 
deformed with four inward electrons forming a square close to the border. 

The pursuit of the FS can also be viewed by observing the compo- 
sition of the wave functions. For the (81)i , the BF with the largest weight 
(35.4%) is the {0-9,11-13}, in which the electrons are divided into two com- 
pact bunchs, and therefore supports directly the (10) c -3 structure . For the 
(82)i , the BF with the largest weight (33.4%) is the {0-8,10-13}. 

When L > 84, the ~N out of the FS is > 6. Since the outward electrons 
would separate (roughly by r u ) from each other, a larger N out definitely leads 
to a larger ring. Consequently, the outward electrons are no more close 
to the core, and the ring-core structure becomes explicit. This is shown in 
Fig. 3c to 3f for the (84) 1 to (90) 1 states, where the outward peak becomes 
larger and larger. The point a separating the inner and outer regions can 
be well defined. Accordingly, the quantities related to eq.(6) to (8) can 
be calculated as listed in Table 2. In particular, a quantity related to the 
uniformity of rotation 

= W/^hl = u out /uj m (17) 

J-out J-xri 

is defined and is also listed. 

It is exhibited in Table 2 that, in the range 84< L < 90, a and iVj n 
are decreasing . This coincides with the reduction of the core of the FS. In 



13 



particular, the Nj n of the FS are one-to-one close to the iVj n from calculation. 
This fact confirms that the FSs are pursued by the first-states. In general 
the N in and l in deviate more or less from those of the FS, this is due to the 
mixing of the FS together with other minor structures ( the inner electrons 
may occasionally go out ,or the core may get slightly excited). E.g., the wave 
function of the (87) i has N in = 3.85 and l in = 7.21 , while its FS has N in =4 
and lin—(hn)b =6 (incidentally, a core-excitation may cause a big increase of 
li n )• Furthermore, the 7 are close to the unity, it implies that the rotation is 
roughly uniform. However, the slight deviation of 7 implies that the system 
is not entirely rigid. 

It is recalled that the P2 of the L < 83 first-states appear as liquid- 
like. However, when N out is neither very small nor very large (say, 6 <N out < 
10 ), the U of the core-ring structure is lower, and thereby the associated FS 
becomes more dominant. This would lead to a clear crystal-like picture as 
shown in Fig.4c to 4g, where the outward electrons form a regular polygon. 
The number of vertexes (from 6 to 10) is just equal to the N out of the FS. 
This fact once again demonstrates the pursuit of the FSs. In general, the 
crystal-like structure can be seen more clearly if p 3 is observed as shown in 
Fig. 5c. 

When N out is larger than 10, due to the rapid increase of U, the asso- 
ciated (Nj„)— N out structure is no more dominant, and therefore the crystal- 
like picture becomes ambiguous again due to the mixing of structures. This 
is shown in Fig.4h and 4i. 

6, ELECTRONIC STRUCTURES (91 < L < 101) 

Inserting f = 1 into eq.(16) and using the same argument as before 
, it is straight forward to prove that the CFS of the L = jV( -^~ 1 - > + N = 91 
states include all the hollow (N in ) h — N out structures ranging from N in = 
to 12 . Therefore the (91) 1 would be nodeless if the core is hollow. On the 
other hand, if the core is inert, all the (N iri ) c -N out are not the CFS (except 
the (12) c -l). Thus, the (91)i is expected to be hollow . This suggestion is 
confirmed by Fig.3g. Similar to the (78)i, the (91)i is also mainly contributed 
by a single BF {1-13} with the weight 82.0% . Due to the lack of coherent 
mixing, the (91)i is liquid-like as shown in Fig.4j. 

For the first-states with 92< L < 101 , we have 

(i) The core-ring structure is explicit as representatively shown by 
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the p\ plotted in Fig.3h to 3j. However, the core may be excited and the 
probability of an electron staying at the center is smaller (3i and 3j). 

(ii) It is noted that a state with a large L would pursue a larger 
moment of inertia to reduce the rotation energy. Since the structures with 
N ou t< 7 have a smaller moment of inertia, these structures are never found in 
the first-states with L > 93 . Specifically, the (92)! is found to have N out =7 
as shown in Fig.4k. 

(iii) Each of the (92) i , (94) i , (96) i , and (98) i states has only one 
FS, this FS has an appropriate N out , and has 9 /G opt ~ 1. Therefore 
these FSs are competitive and are expected to be dominant. This point is 
confirmed by the associated p-i (cf. Fig. 4), where a crystal-like picture with 
the N otii — side polygons is seen. Furthermore, the N in of the FS of the above 
four states are 6, 5, 4, and 3 (cf. Table 2), while the N in are 6.26, 5.01, 
3.89, and 2.91, respectively . These values are one-to-one close to each other. 
Thus, the pursuit of the FSs is further confirmed. Besides, the FS of the 
above four states have U n = (k n )b, namely 15, 10, 6, and 3 (cf. Table 2), 
respectively. The corresponding l in calculated from p\ are 19.70, 12.51, 7.38, 
and 4.19, respectively. The latter set are always one-to-one bigger than the 
former set due to having a slight core-excitation. 

(iv) When L > 100 , the excited core (i.e., > (li n )b ) begin to 
compete seriously with the inert core. For the L = 100 states, the competing 
FSs are the (4)^-9 and (2) c -ll as shown in Table 2. The 9 /G opt of the former 
(latter) is considerably smaller (larger) than one. It is noted that, when L 
is large, the outer ring would shift a little outward to increase the moment 
of inertia to reduce the rotation energy. Thus, a small increase of 9 is of 
advantageous, while a decrease of 9 is not. In fact, it is the (2) c -ll wins in 
the competition and is pursued by the first-state, while the (4)^-9 is pursued 
by the second-state. This is shown in Fig.3h and 3i, and in Fig.4o and 4p. 
For the (101)i, the (3)^-10 is the only FS, and is expected to be dominant 
as shown in Fig.3j and 5d. 

(v) For the (93) i , (95) i , and (97) i , the N out of their FS are smaller 
than 7 and therefore is not competitive. Although the (99) i has N out =8, 
however its 9 /G opt is quite small. Thus these four states do not have a 
competitive FS, and therefore do not have a clear-cut geometric structure to 
pursue. They are liquid-like. Nonetheless, their p\ are more or less similar 
to Fig.3h, thus they still have clear core-ring structures. 

(vi) All the first-states with 92< L < 101 rotate uniformly, they have 
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7 ~ 1 except the (99)i and (101)i, The FSs of these two states have a con- 
siderably smaller g /G opt (cf. Table 2). Thus, due to eq.(14), if they rotate 
uniformly the ring would be too close to the core . To avoid being too close, 
u out would decrease a little. In this way, although the rotation energy may in- 
crease a little, the potential energy may thereby considerably decrease. This 
suggestion is confirmed by the fact that their 7 is really smaller. Inciden- 
tally, since the angular momentum l out is strongly constrained by symmetry 
via eq.(9a), and therefore can not be adjusted freely , the decrease of uj out 
would cause an increase of I out via the relation l out = I ou t0J O ut- 

7, ELECTRONIC STRUCTURES (L « 200) 

The main finding of the above study is the pursuit of the FSs . Does 
this experience work when L is much larger? To clarify this point we shall 
no more go to the states one-by-one, instead we choose arbitrary a range 
196< L < 201 for the studying. Let us first evaluate the accuracy of the 
calculation in this range. E.g., the energies of the (199)i state calculated 
with 6000, 9000, and 12000 BFs, respectively, together with the a, jVj„, U n 
and 7 are listed in Table 3. One can see that, although the convergency is 
not very good, it is qualitatively acceptable. 

Table 3 The energies and the quantities extracted from the p\ of the 
(199)i in accord with the increase of the number of BFs. 



Number of BFs 


a 


N in 




7 


£((199)0 


6000 


3.121 


3.96 


10.56 


1.268 


739.97 


9000 


3.123 


3.97 


10.75 


1.238 


739.86 


12000 


3.125 


3.98 


10.86 


1.227 


739.81 



The FSs are shown in Table 4. The FSs with the core inert (k = 0) 
are found to have a too large g/ G opt , and therefore are not listed. Whereas 
an excited core is pursued. On the other hand, a highly excited core (l\ > 6) 
would lead to a too small g/G opt as shown in the table. Thus, too weak and 
too strong core-excitation are both not appropriate. 

Table 4 A continuation of Table 2 for the first- states with (196 < 
L < 201). 
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Favorable structures and 




Quantities evaluated from 




their relatec 


features 




the p! of the first-states 


L 


FS 


Zi 


kn 


9 


g/ G opt 


a 


N in 


kn 


7 


196 


(5) h -8 


6 


40 


1.5612 


0.72 










196 


(3) ft -10 


6 


21 


1.5811 


0.59 


3.16 


3.75 


13.88 


1.11 


196 


(2) h -ll 


4 


9 


1.9437 


0.52 










197 


(5) h -8 


3 


25 


2.0736 


0.96 


3.48 


4.80 


24.61 


1.05 


197 


(4) h -9 


5 


26 


1.7097 


0.71 










197 


(3) h -10 


3 


12 


2.1506 


0.76 










198 


(4) fc -9 


3 


18 


2.1082 


0.85 


3.25 


3.90 


18.65 


1.10 


198 


(2) h -n 


5 


11 


1.7581 


0.47 










199 


(4) h -9 


1 


10 


2.8983 


1.20 


3.13 


3.98 


10.86 


1.23 


199 


(5) h -8 


5 


35 


1.7113 


0.80 










200 


(5) h -8 


2 


20 


2.3717 


1.01 


3.47 


4.93 


20.88 


1.10 


200 


(3) h -10 


4 


15 


1.9235 


0.70 










200 


(2) h -ll 


6 


13 


1.6172 


0.43 










201 


(3) h -10 


1 


6 


3.1225 


1.10 


3.00 


3.06 


7.27 


1.24 


201 


(2) h -ll 


1 


3 


3.4641 


0.93 










201 


(4)^-9 


6 


30 


1.5916 


0.66 











For the (196)i none of the FSs are superior (their g /G opt are too 
small), therefore this state would have a strong mixing of structures and 
would be liquid-like. Among the three FSs of the (197)i, the (5)^ — 8 has 
its g closer to G op t , thus this FS is predicted to be dominant. Similarly, 
based on Table 4 , the (4)^—9 is predicted to be dominant in (198)i and 
(199)i ,the (5)^—8 is predicted to be dominant in (200)i, and the (3)^—10 
is predicted to be dominant in (201)i. It turns out that, for the case with 
a dominant FS, the predictions are nicely confirmed by the calculation. E.g., 
the N in of the above FSs of the (197)i to (201)i are 5, 4, 4, 5, and 3, while 
the corresponding N in extracted from p 1 are 4.80, 3.90, 3.98, 4.93, and 3.06 . 
The l in of the above FSs are 25, 18, 10, 20, and 6, while the corresponding 
l in extracted from p x are 24.61, 18.65, 10.86, 20.88, and 7.27 . These values 
are amazingly one-to-one close to each other, and thus the analysis based on 
the FSs is convincing. Furthermore, the associated p-i and confirm also 
the predictions. Representative examples are given in Fig. 4q, 4r, 5e, and 
5f. 

It is noted that the (199)i and (201)i have a considerably larger 7. 
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On the other hand, their most competitive FSs have a larger g /G op t. Thus, 
if these states rotate uniformly, the ring would be too far away from the core 
(cf. eq.(15)). To avoid being too far away, the ring rotates a little faster to 
reduce the moment of inertia without altering l out . This is the reason why 
they have a considerably larger 7. 

In general, when L is large, the size of the system would increase, 
the core-ring structures become more clear-cut. Besides, the core would have 
a higher excitation. As a result, all these states are hollow as shown in Fig. 
3k, 3Z, 5e and 5f. 

8, MAGIC NUMBERS 

The above discussions demonstrate that, based on the FSs, the 
structures of the first-states can be more or less predicted. In this section we 
shall see that the energies are also strongly related to the FSs. Let us go back 
to Fig. 2 where platforms and shoulders are shown. A platform starting at L a 
and ending at Lb = L a + 1 implies E((Lb)i) = E((L a )i)+tiw , i.e., the (£&)i is 
an cm. excited state of the (X a )i- This fact implies that the internal energy ( 
the energy without the cm. motion) of the (L b )i states are relatively higher. 
This is also the case if a shoulder appears. In this case, L a is a candidate 
of a magic number (CMN). Evidently, if the (L a ) 1 has a competitive FS and 
the {Lb)i does not have, a CMN arises. For example, the (78) 1 is inherently 
nodeless and is able to get access to all symmetric configurations , while the 
(79)i has only one CFS (12) c -l which is unfavorable to binding. Thus the 78 
appears as a CMN. Similarly, the (91) 1 is inherently nodeless (if the core is 
excited), while the (92) 1 has only the (6) c -7 ( which is not competitive due 
to N out = 7), thus 91 is a CMN. The (lll)i has a competitive FS (3)^-10 . 
Although the (1 12)i has two FSs, namely the (6)^—7 and (5)^-8 , however 
the former has a small N oui while the latter has a too small 9 /G opt = 0.78 . 
They are both not competitive, thus 111 is a CMN. Finally, The (1 18)i has 
a number of competitive FSs, namely the (3) c -10 , (5) c -8 , and (4)^-9 , while 
the (119)i has only one FS (6)^-7 , which is not competitive due to having 
N out = 7. Thus 118 is a CMN. These examples exhibit that the CMN can be 
more or less predicted. 

9, SUMMARY 

The electronic structures of the first-states have been studied. By 
an analysis of symmetry constraint and by performing numerical calculation, 
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we have obtained a clear picture of the core-ring structures. When L is small 
(78< L < 83), the core and ring are connected. When L is larger than 83, 
the core-ring structure becomes more and more explicit. When L < 100, 
the core remains inert (the (91) i is an exception). When L is larger, core 
excitation begins to compete. When L is much larger (say, L m 200), core 
excitation becomes dominant and the states are hollow. 

The number of particles and the amount of angular momentum con- 
tained in the core (ring) are not only determined by dynamics, but depend 
seriously on symmetry constraint. Due to the constraint, for a given state, 
it is advantageous to pursue a specific kind of structure, but disadvantageous 
to pursue another kind. This leads to Table 2 and 4, where the favorable 
structures (FSs) of each state are listed. 

The identification of the FSs is the main result of this paper. Based 
on the FSs, the structures of the first-states can be predicted to a great 
extent, the formation of crystal-like structure and the appearance of magic 
numbers can be explained. In particular, if a first-state has a remarkably 
competitive FS (both the N out and g /G opt are appropriate), the FS would 
be pursued, and the state would be crystal-like and possess the associated 
(Nin)— N out structure. If the L = L a states contain one or more than one 
competitive FSs while the L = L a + 1 states do not contain, then L a is a 
CMN. 

The FSs can provide us an objective base for the further classification 
of states. The states having the same FSs can be grouped into a kind, e.g., 
all the L = 87 , 96, 105, • • • contain a single FS (4) c — 9, thus they belong to 
the same kind and their first-states would have the same (4) c — 9 structure. 

Although only a N=13 system is concerned in this paper, the idea, 
the way of analysis, the qualitative results are quite common to the 2- 
dimensional systems with an attractive center. In fact, both this paper and 
the previous ref.[8] provide qualitatively similar message. Thus, it is not 
doubted that the physical picture provided by these two papers can be gen- 
eralized to the systems with an even larger N. Where, the identification of 
the FSs is again a key to understand the electronic structures . 
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Fig.l pi(r) of the first-state (82)iwith 6000 (a), 9000 (b), and 12000 (c) 
basis functions. The unit of length in this paper is Jh/m*LU = 194.71 A ■ 



Fig. 2 e(L) as a function of L. frc<Jo=3meV is assumed. 
Fig. 3 pi(r) of the first-states (Fig.3i is for a second-state). 
Fig. 4 The contour plot of the two-body densities p 2 (r, r 2 ) as a function 
of r. The given r 2 is marked by a white spot. The lighter region has a larger 

P2 ■ 

Fig. 5 The contour plot of the three-body densities p 3 (r, r 2 , r 3 ) as a func- 
tion of r. The given r 2 and r 3 are marked by two white spots. Refer to 
Fig.4. 
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